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I. INTRODUCTION 

In the study of systems of many interacting particles, one of two situations might arise. If the system is isolated, 
or if the interaction with its environment allows it (is invariant under the reversing of time, for instance), one could 
observe an equilibrium state, in which the probability of any configuration is simply related to the energy of that 
configuration via the Gibbs-Boltzmann law. In principle, from that information, and the expression of those energies, 
one could calculate any equilibrium quantity that might take his fancy. Obviously, those calculations could still be 
extremely difficult, but the general framework provided by the Gibbs-Boltzmann law systematically solves the first 
part of the problem. In the other situation, where no equilibrium state can be observed, there is no such framework. 
The system might reach a non-equilibrium steady state, where the probabilities of the configurations do not depend on 
time, but there is no way to calculate those probabilities a priori, and one has to solve the whole dynamical equation 
that describes the evolution of the system to obtain the desired quantities. 

There have been, however, many attempts to generalise the concept of Boltzmann weight to non-equilibrium 
systems. The best candidates to play the role of the energy in these systems are large deviation functions, which 
contain information on the probabilities of rare events or configurations, in the limit of some large extensive quantity 
(usually time or system size) [21 |4H9] . Those large deviation functions have for instance helped uncover some very 
general symmetries, called 'fiuctuation theorems', that are verified by systems however far from equilibrium |121 113| . 
The study of large deviation functions is therefore an important task to statistical physicists. 

Many of those systems that display non-equilibrium steady states describe the transport of carriers (e.g. of mass, 
electrical charge, or thermal energy), that may interact with each-other, through some domain (a one-dimensional 
channel, for example), and driven by an external field in the bulk of the system, and/or unbalanced reservoirs 
at its boundaries. One may think, for instance, of a metallic wire conducting electrons between two masses at 
different potentials, or an artery conducting blood cells between two organs at different pressures. Because of those 
driving forces, the system exhibits a non-vanishing macroscopic current in its steady state. That current is related 
to the microscopic entropy production that comes from the breaking of detailed balance and time reversal invariance, 
and which is characteristic of non-equilibrium states (in some cases, that relation between macroscopic current and 
microscopic entropy production is a very simple one, as one can see in [12' and appendix [A|). 

One of the simplest examples of driven particle models, and one of the most extensively studied, is the Asymmetric 
Simple Exclusion Process (ASEP). It consists of a one-dimensional lattice, the sites of which hold particles that jump 
forwards and backwards stochastically. The particles interact via hard-core repulsion, so that there can be only one 
particle on a given site at a given time (hence 'exclusion'). They jump preferentially to one direction, which accounts 
for the driving force in the bulk of the system (which makes it 'asymmetric'). Each side of the system is connected to 
a particle reservoir characterised by a fixed density. The ASEP has many qualities which explain the extent to which 
it has been studied in the past twenty years or so [2 [71 HI IHl IT^HTOl 1511 [55] . First of all, it is simple in its definition, 
and has the mathematical property of being integrable, which makes it a good candidate for analytical calculations 
and exact solutions. This integrability property implies that the methods used to treat the model are usually not 
general and transposable (except to other integrable systems), but the actual results could give us valuable insights 
into the general behaviour of generic non-equilibrium systems. Moreover, the ASEP has connexions with many and 
various other systems, such as ribosomes travelling on a m-RNA strand [H [^nk22j (which is what the ASEP was 
originally meant to describe), random polymer models [23], growing interfaces [10ll2l], pedestrian and car traffic [52) . 
quantum spin chains [JS] , random matrices p5H?f] and even pure combinatorics [IHl [SD] • 

In the present paper, we add our own effort to the long history of results on the steady state of the ASEP and 
the fluctuations of the current it exhibits, by solving the long-standing problem of obtaining the distribution of those 
fluctuations in the open case. Because the system is out of equilibrium, the choice of the boundary conditions matters 
greatly, much as for systems with long-range interactions. The statistical ensembles are not equivalent, and fixing the 
number of particles (as in the periodic case) or not (as in the open case) makes the results, as well as the methods 
applicable to their acquisition, significantly different. For instance, the full current fluctuations for the periodic 
TASEP (where the T stands for 'totally', i.e. the particles jump in only one direction) were obtained by Bethe Ansatz 
in [35| . and the same method was used many years later for the periodic ASEP in |36| . but that method turned out 
to be inapplicable to the open case, precisely because of the non-conservation of the total number of particles. The 
Bethe Ansatz could also have led to the distribution of the steady state probabilities in the periodic case, were it 
not trivially flat, but the same could not be done in the open case, and another method had to be devised, namely 
the 'matrix Ansatz' pi I33j. by which some recursion relations 32J that were found in the weights of conflgurations 
of systems of different sizes are encoded in algebraic form. A flrst extension of that method was used to obtain the 
diffusion constant (second cumulant of the current) in the open ASEP [3j. During the last couple of years, the author 
and collaborators managed to generalise that method further to calculate all the cumulants in the open case, flrst for 
the TASEP |32] , then for the ASEP [33] . In both cases, the results are exact for any values of the parameters and for 
any finite size of the system. However, at the time, large portions of the proof were guessed rather than fully worked 
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out, so that the results were given as a conjecture along with numerical evidence to support them. In this paper, we 
give the complete algebraic proof of the validity of our Ansatz, and explain how it gives us access not only to the 
fluctuations of the current, but also to any spatial observable conditioned on the mean current. We also show how 
the Ansatz extends to the periodic case, and to the spin-| XXZ chain with non-diagonal boundary terms. 

The layout of the paper is as follows: in section [llj we define the model, and do a quick review of some known 
steady state properties; in section |III[ we restate the problem of finding the fluctuations of the current as that of 
finding the first eigenvalue of a deformed Markov matrix, and we define the 's-ensemble' as the distribution of the 
principal eigenstate of that same matrix; in section [TV] we present our perturbative matrix Ansatz, along with its 
proof (the technical portions of which are carried out in the appendixes) , an alternative formulation which was used 
in [43], and its equivalent for the periodic case and the XXZ chain; finally, in section |vj we give a quick overview of 
how the explicit calculations of the cumulants of the current were carried out using our Ansatz. 



II. DEFINITION OF THE MODEL AND STEADY STATE PROPERTIES 



A. The open ASEP 

The open Asymmetric Simple Exclusion Process (ASEP) in continuous time is defined as follows. Let us consider 
a chain of L sites, numbered 1 through L, each of which can hold at most one particle. The occupation of site i is 
denoted (equal to for an empty site, and 1 for an occupied site). Those particles can jump stochastically to the 
right with rate p = 1 and to the left with rate 9 < 1, provided that the receiving site be empty. In addition, each end 
of the chain is connected to a reservoir of particles, so that particles may enter the system at site 1 with rate a or at 
site L with rate 5, and exit the system from site 1 with rate 7 and from site L with rate j3 (fig. [l]) (the rate p can be 
set to 1 without any loss of generality). 

a 1 /3 
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FIG. 1. (colour online) Dynamical rules for the ASEP with open boundaries. The rate of forward jumps has been normalised 
to 1. Backward jumps occur with rate g < 1. All other parameters are arbitrary. The jumps shown in green are allowed by 
the exclusion constraint. Those shown in red are forbidden. 



At any time i, the system can be described by the probability vector ||Pt)), of which the entries Pt{C) give the 
probability of being in the configuration C = {Ti)i..L at time t. This probability depends only on the initial condition 
||Po))i and verifies the master equation: 

^ IIP,)) = M||P,)) (1) 



solved by 

where the Markow matrix M is given by: 



dV 

||P,))=e*^'^||Po» (2) 



L-l 



M = Mo + ^ M, + Mi (3) 

1=1 

with Mi containing the jumping rates between sites i and i + 1, and Mq and Ml corresponding to the couplings with 



the left and right reservoirs (see equ.(14| with /i set to for the explicit expression of those matrices). 

Each non-diagonal entry of M contains one of the aforementioned rates, and is non-zero only if the initial and final 
configurations differ by the jumping of exactly one particle. The diagonal entries contain (minus) the rate at which 
the system leaves a given configuration, which is also the sum of the rates from that configuration to any other, so 
that the sum of each column of M is (i.e. the evolution of the system conserves probability, and M is a stochastic 
matrix). From the Perron-Frobenius theorem, we know that the largest eigenvalue of M is 0, and that, for large times, 
||Pt)) converges to the corresponding right eigenvector ||P*)) (the corresponding left eigenvector is the vector with all 
entries equal to 1, and will be denoted by ((1||). 
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B. Matrix Ansatz for the steady state and phase diagram 



The exact expression of the steady state can be written in the form of a matrix product state (also caUed 

'matrix Ansatz') [33j : 



P*(C) = ^{W\fl inD + (1 - n)E) \V) , (4) 

where the matrices D and E, and the vectors {W\ and \V), are defined by the following algebraic relations: 

DE-qED^{l-q){D + E) 
{W\{aE~^D) ^ {l-q){W\ 

{PD-6E)\V)^{l-q)\V) . (5) 

and Zl = {W\{D + E)^\V). 

Let us remark here that throughout this paper, we use the standard bra/ket notation for vectors from the space on 
which matrices D and E act (e.g. \V)), and the doubled bra/ket notation for vectors from the configurations space 
(e.g. IIP*))). 

As stated in (Hi), the weight of a given configuration C = (T,i)i in ||P*)) can be written as the product of L matrices 
D or i?, contracted between two vectors. The ith matrix in the product corresponds to the occupation of the ith site 
in C: it is D if Ti = 1 and _E if = 0. The first algebraic relation in ([s]) encodes a recursion between the steady state 
probabilities of systems of successive sizes. Combined with the two other relations, it allows, in principle, to compute 
explicit expressions of any of those probabilities, although that computation would be extremely impractical. However, 
it can be used, in a most elegant manner, to compute the mean values of certain observables, like the particle current 
passing through the system, or the local density ;32j, which are more physically relevant than the probability of a 
single configuration. For instance, for a system of size L, the stationary current J is given by J = {l~q)Zi^_i/Zi^, and 
can be expressed in terms of certain orthogonal polynomials [23 EE] ■ Those calculations are especially easy to carry 
out in the simpler case of the Totally Asymmetric Simple Exclusion Process [S^ (TASEP, for which g = 7 = (5 = 0) , 
and even more so in the case a = /3 = 1 [30 . 

In the large size limit L — cx), those quantities (mean current J and mean density p{x) where x = i/L) can be 
used to characterise the different phases in which the system may find itself, depending on the values of the boundary 
parameters. Those phases are best described in terms of the effective reservoir densities 

p, = l/(l + a+) (6) 
p, = b+/{l + b+) (7) 

where 



(1 - (7 - a + 7) ± ^(1 - (? - a + 7)2 + 4a7 
a± = ' (8) 

{l-q-p + 5)±^{\-q- P + 

= ■ 

We may note that those densities verify the relations aj^—^j^zTf^ = (l^?) ^ind (izr^^*^^ = (1^9)i which can 
be related to the second and third equations in ([5|. 

The phase diagram of the open ASEP [H [14] , which is essentially identical to that of the TASEP up to a factor 
(1 — q) in the current, is divided in three regions (fig. [2]): 

• Maximal Current phase (MC): pa>^/'2. and pb < 1/2 (i.e. a+ < 1 and 6+ < 1) ; the mean current is J = ^^^^^ ; 
the mean density is 1/2 in the bulk, with an algebraic crossover to the boundary density at both ends. 

• Low Density phase (LD): pa < 1/2 and pb < I — Pa (i-c a+ > 1 and 6+ < a-|.) ; the mean current is J = 
(1 — q)pa{^ — Pa) ; the mean density is pa in the bulk and at the left boundary, with an exponential crossover 
to Pb at the right boundary. 

• High Density phase (HD): pa > I — Pb and pb > 1/2 (i.e. 6+ > 1 and a+ < 6+) ; the mean current is J = 
(1 — q)pb{^ — Pb) ; the mean density is pb in the bulk and at the right boundary, with an exponential crossover 
to pa at the left boundary. 
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FIG. 2. (colour online) Phase diagram of the open ASEP in terms of the boundary densities pa and pb. The profiles in each 
phase show the schematic behaviour of the mean density in the system for large sizes. 



In order to access the fluctuations of the current, and not only its mean value, we need to find the steady state 
probabilities as a function of the time-integrated current as well as of the configuration. That is what we propose to 
do in this paper. 



III. CURRENT-COUNTING MARKOV MATRIX AND THE S-ENSEMBLE: 



Suppose that we want to keep track of the number of particles that jump over the bond that links sites i and i + l 
through the evolution of the system, starting from some initial state l^b))- One easy way to do this is to multiply the 
off-diagonal entries of Mi by a fugacity e*^' , so that every time those jumping rates are used in the evolution of the 
system, the weight of the corresponding history gets multiplied by e^' if the jump was made forwards, or e"''' if it 
was made backwards. After a time the weight of any history C{t) carries an extra weight e"^*^'''-*-'-'''' , where J{C(t)) is 
the (algebraic) number of particles that jumped from site i to site which is precisely the time-integrated current 
that went through that bond. One can then access the moments of that current simply by taking derivatives with 
respect to pii. 

In general, one can do the same with any and all of the bonds, using different fugacities. The corresponding Markov 
matrix can be written as 22 [SU] : 



L-l 



(10) 



where 



-a 76 



■0 





0" 


-q 


e^' 





ge-^' 


-1 





.0 









-S |3e^'^ 



(11) 



(it is implied that Mq acts as written on site in the basis {0, 1} and as the identity on the other sites, and the same 
goes for M]^ on site L; similarly. Mi is expressed by its action on sites i and i+1 in the basis {00, 01, 10, 11} and acts 
as the identity on the rest of the system). 
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By using this deformed Markov matrix in the time evolution of \\Ptj) instead of the usual one, one gets 

||P,({/.,})))=e**^^H>||Po)) (12) 

which is the Laplace transform of the joint probabilities of the configurations and the time-integrated currents, with 
respect to these currents, at time t. Taking a fcth order derivative in any of the /ijS and then projecting this vector 
onto ((1|| yields the fcth moment of the corresponding current, up to a normalisation. 

In the long time limit, the matrix e**^*^*' converges to the projector onto its principal eigenvector \\P\fii})), with the 
eigenvalue e*^^^*^-* where is the largest eigenvalue of M^.j, so that 



pts({^.}) 



(13) 



and E{j_ii} is therefore identified as the generating function of the cumulants of the instantaneous currents Ji/t, which 
is the Laplace transform of the large deviation function of those same currents, as explained by the Donsker-Varadhan 
theory of temporal large deviations [2 [SD]- That is the main quantity that we want to calculate. Our problem thus 
reduces to that of finding the largest eigenvalue of the current-counting matrix M^.j . The corresponding eigenvector 
II P{^.})) also holds important information, as we will see below. 

Let us first make things a little simpler by noticing that one can go from any set {fij} to any other set {fi'^ by a 
matrix similarity, as long as X^i^o ~^f=o (^^^ Appendix A for a detailed derivation of that statement, and 

[T2j). That means that the eigenvalues of Mj^.j only depend on /i, regardless of how the fugacities are distributed: 
the currents through each of the bonds are all exactly equivalent. In particular, there is a set for which fi 

^ where A is the quantity conjugate to the entropy production in the system (see appendix 



A log 



A). This is 



an easy way to prove the Gallavotti-Cohen symmetry for the current in this system, and it shows that fne current 
and its fiuctuations are simply proportional to the entropy production. That entropy production being non-zero is 
the defining characteristic of a non-equilibrium system. 

All this being said, we can now consider, without any loss of generality, the case where only the first bond (between 
the leftmost reservoir and the first site) is marked: /zq = , fJ-i^o = 0, so that the individual jump matrices we will 
work with are given by 



Mo(Ai) = 



-a 7e 



M,, = 



0" 

1 

g -1 





, Ml = 



-6 /?■ 
6 -/3 



(14) 



and that 



with 



IIP I 



e*^('^)||P„ 



(15) 



fc=i 



the exponential generating series of the cumulants of the current. 
The large deviation function for the instantaneous current, G{j) 



i log[P(J/t = j)], is related to E{^i) by 



(16) 



G{j) = fij - E{^^) 



(17) 



Moreover, one can show [39] that the principal right eigenvector ||P;^)) of M^, and its principal left eigenvector((Pyj||, 
hold (respectively) the probabilities of observing configurations coming from the steady state, and the probabilities 
of having come from configurations while in the steady state, all conditioned on having observed a mean current 



P,iC)=P{{C,t} {P\t = -^} & J^^Eif,) 
P^{C)=P{{C,t} {P\t = +^} & J = -^E{f,) 



(18) 
(19) 
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so that the product of the two is the probabihty of observing a configuration at any time, conditioned on the mean 
current (up to a normaUsation) : 

P,(C)P,{C)=p(c\j = ^E{f,)) (20) 

The ensemble defined by these probabihties, with /i as a parameter, is sometimes called the 's-ensemble' |^ (the 
reason being that fj, is often noted s), and contains all the information needed to build the joint large deviation 
functions of the current and any spatial observables (i.e. depending only on C). We will now construct a matrix 
ansatz that holds the exact expressions of those probabilities and of E{^), as series in /x, up to arbitrary orders. 

IV. PERTURBATIVE MATRIX ANSATZ FOR THE S-ENSEMBLE: 

In this section, we will show that we can define two transfer matrices and Uf^, such that: 

[M^,C/^T^]-0 (21) 

ro-[l]c.C'-||l))((l|| (22) 

C/o||l)) = IP"-)) (23) 
where the weights of and are expressed as products of matrices between two vectors, much as in (|4|. Moreover, 



equ. (21) for /i = allows to recover that original Matrix Ansatz, as we will show below. 

Those three relations will be used to prove that the transfer matrix UpT^ is almost a projector onto the leading 
eigenstate of M^, and that when applied repeatedly, the precision in orders of /i of the projection increases. In other 
terms: 

{U,T,f ^\\P,)){{P,\\+0{^i'') (24) 

up to a multiplicative constant of order 1. 

We will then use that relation to show the main results of this paper, namely that 

\\P,))^^){UpT,nP*))+0{^,'■+^) (25) 



((^mII = -^^m{u,T,Y + o{^,''+') (26) 

where ^ ^ {{l\{UpTpf\P*)) , and that 

_ {{l\\M,{UpT,nP*)) 

mm,T,np*)) )■ 

Those results hold for any integer fc, so that, in essence, we have complete exact expressions for {{Pfi\\ and 

E{^i)^ expanded as infinite series in /i. 

A. Definitions and commutation relations: 

Let us consider two matrices d = D — 1 and e = E — 1 where D and E are defined as in ([s]), and a matrix A^, such 
that: 



de — q ed— {\ — q) 

Apd = e'^ dA^ . (28) 
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The first of these relations, which is a simple consequence of ([5|, defines the algebra of a q-deformed harmonic 
oscillator j54j . of which e is the creation operator, and d the annihilation operator. 
Let us also define two more boundary vectors {W\ and \V) by: 



m - rf) - 6(1 - e)] \V) = 
(W^|[a(l-e)-7(l-d)] =0 



and let us recall that 



[P{l + d)-S{l + e)]\V) = il-q)\V) 
{W\[a{l + e) - 7(1 + d}] ^ (1 - q){W\ 

By writing 

-''^o^o = ^1.1 = 1 ; -^^1,0 — d , Xq^i — e 
we can finally define the weights of and between configurations C = (TjQl.l and C = {Ti)i,,L' 



1 ^ 

XC,C')^—{W\A,Y[X,^,^,\V) 



(29) 



(30) 



(31) 



(32) 



i=l 



with Zl = {W\{2 + d + e)^\V), and 



T^{C,C')^{W\A^l[Xr^^,,\V}. 



(33) 



Those weight are entirely determined by the algebra defined above. Specifically, one can obtain the weights of 
buy using ( 28 ) and ( 30 1 , and those of buy using ( 28 ) and ( 29 1 . 

We may note that the matrix Afj_ is set between the left boundary vector and the first matrix because it is the bond 
between the left reservoir and the first site that is marked. For a general set of weights {fii}, we would have to add 
matrices A^^. between X^.y and ^Ti+i.T^/^i in both of the products above, so that: 



1 ^ 
;7{^,}(C,C') = -^{W\A,,l[Xr^yA^AV) 

L 

T^^,}{C,C') = {W\A,,l[Xr^yA^^\V). 



(34) 
(35) 



In appendix [Bj we derive equation (21 1, using a method closely related to the matrix Ansatz for niultispecies ASEP 
on a ring [37j, and which makes use of the so-called 'hat matrices'. That equation is the main point to our Ansatz: 
the transfer matrix U^T^ that we built has the same eigenvectors as M^, so that we can try to extract the information 



we need from it instead of M,, 



For /i = 0, one particular solution to (28 1 and (29 1 is d = e = = 1 
To(C,C') = which we can set to 1. This proves (22 1. Furthermore, projecting Uq onto 

over all configurations C in ( 32 1 , so that 



so that for any C and C, we have 
1)) means summing 



iC\\Mn=Y.^,{C,C') ^ (Ty|Aon(^n.o + ^na)|^) 



(36) 



We can set Aq to \, and remark that for = 0, we have {X^fi + Xn^i) = 1 + e = E and that for = 1, we have 
{Xnfi + ^Ti,i) = d + 1 = I?, SO that this expression is exactly that of P*{C) as given in (j4|, which proves (23). 
Using relations (21), (22) and (23) together, at fi — 0, we get 



[MM] = = (m||p*))((i||) - (||P*))((l||Af 



(37) 
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Since we know that = (because M is a stochastic matrix), this impHes that M\\P*)) = 0, which yields the 

original matrix Ansatz Q [33] . 

This alternative proof of Q relies on the fact that the transfer matrix C/^T^ is a projector in the limit /i — > 0. 
It would be interesting to determine whether for other situations with matrix product states, one can generically 
find a transfer matrix that commutes with a deformation of the dynamics of the system and is a projector in the 
non-deformed limit. One could for instance look at the ASEP in discrete time with different versions of the update 
[3H] , or at the multispecies ASEP on a ring [37] . 



Validation of the perturbative matrix Ansatz: 



To prove (251 and (27), we use the relations derived above. Since, for /i = 0, the matrix UqTo is the projector onto 
the principa 



eigenspace of M, one can write, for an infinitesimal /i: 



(38) 



where A^ ~ 1 + O (/x) is the largest eigenvalue of U^T^^, and '-^ O {p) is the part of U^T^ that is orthogonal to its 
principal eigenspace, and has eigenvalues of order ^. In other words, UfjT^ is almost a projector, with an error of 
order /i. 

Since r^\P^)) = and = 0, one has that: 



{U,T,f = Kl\\P,)){{P,\\+rl (39) 

so that the difference from the projector onto \P^^'f|{{P^l\ is now ^ O {iJ.'^)- 

Let us now remark that the parts of \\P*)) and of ((1|| which are not in the kernel of are of order fj,, so that both 
r^\\P*)) and {{iWr'^ are of order /i'^"'"^. It follows that {U^Tf^)''\\P*)) is proportional to with an error of order 
(and the sa me g oes for ((/^||), which proves (25). 



Equation (27) is then proven by simply applying Af^ to (C/^T^) 



{l\\M,iU,T,)''\\P*)) = Ei^,)Al{{l\\P, 



(l||(C/^T^)^m)=A^((l||P^ 



ip,\\p*l 



(i||M^'^;:iip*) 



ip,\\p'i 



(i||<ii^*> 



(40) 
(41) 



where imM^r'^jP*)) is of order ^'^^^ because is of order n and r^||P*)) is of order fi'^'^^ , and ((l||r^||P*)) is of 

order /x'^^^ for the reason given above. The ratio of those two equations is therefore equal to E{^) up to order 



C. Formulation as a matrix product: 

The formulation we gave here of the perturbative matrix Ansatz in terms of transfer matrices is quite different from 
that which was given in [33]. They are of course equivalent, which is what we will show in this section. 

The main point that needs to be made here is that, unlike the Markov matrix, which is a sum of elementary 
matrices, the transfer matrices Up, and are products of the elementary matrices Xt-.^j^t-.,^^ , so that a product of 
those transfer matrices can be seen as a tensor network, the tensors being of order four: (X^.^^ ,-.,^^)^ where i and 

j are the internal indices of X (fig. |3|. 

A consequence of this is that the object (C/^T^)'=||P*)) can be written in terms of the columns of this tensor network 
instead of the lines. Let us therefore define, by recursion (and denote the product between successive rows by a tensor 
product (g)): 

Dk+i = {l(Sl + d<^e)(^Dk + {l(^d + d(^l)<^Ek 
Ek+i = (1 «) 1 + e (g) d) (g) Pfc + (e (g) 1 + 1 (g) e) (g) Dfe 

4'=+!) = (g (g 4'=) (42) 
with Dq ^ D, Eo ^ E and 4°^ = 1. and 



\Vk+i) = \V) (g \V) (g \Vk) 
{Wk+i\ = {W\®{W\®{Wk\, 



(43) 
(44) 
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FIG. 3. (colour online) One of the tensor networks that compose Ufj,Tfj,\\P*)) when expanded in terms of the intermediate 
configurations at each step of the product. Lines represent the transfer matrices and T^, whereas columns represent Ek or 
Dk- See the detailed explanation below. 



with l^o) = \V) and {Wo\ {W\. 



In this formalism, equation (25) becomes: 



IIP,, 



{Wk\Al^'> n inDk + (1 - n)Ek) \Vk) + O (m 



k+l\ 



(45) 



i=l 



To give a simple explicit example, in (fig |3| is shown one of the tensor networks which compose U^T^jP*)) for 
L = 5, namely U^{C,C')T^,{C' ,C")P*{C"), with C = {1,0,0,1,1}, C = {1,0,1,1,0} and C" = {0,1,0,0,1}. The blue 
rectangle corresponds to r^(C',C"), the black one to P*{C"), and the red rectangle is one of the elements in Di, 
namely ® (g) E. Summing over the second and third indices in any column gives Ei or Di, depending on 
whether the first (upper) index is or 1. 



While the transfer matrix formulation ( 25 1 is better suited to the algebraic proof of the Ansatz, this matrix product 
formulation (45 1 is useful for doing explicit calculations, like those of the cumulants of the current (see section |v]). Let 
us also note that the matrices Ek and Dk are related to the ones used in the matrix Ansatz solution of the multispecies 
periodic ASEP [37], although for the moment we have no understanding of why that is the case. 



D. Periodic case and XXZ spin chain: 



The same Ansatz can be applied to the periodic case with just one alteration 
products between boundary vectors, one has to take a trace. Moreover 
defined. This can be written as 



instead of projecting the matrix 
only one transfer matrix T'"^'' needs to be 



\c,c') = Ti-[a,1[x,^^,, 



(46) 



if the marked bond is the one between sites L and 1. One can then show that [M^, T^^'"] = and the rest follows, by 
replacing C/^T^ by T^'^'" in every equation, using the steady state IIat)) (with coefficient 1 for all configurations with 
N particles), and making only one tensor product per order in (42). See appendix [C] for all the derivations related to 
the periodic case. 

In appendix p we show how a special choice of the parameters {/i^} allows our Ansatz to be applied to the spin-| 
XXZ chain witn non-diagonal boundaries. We also note that the structure of our construction is strikingly similar to 
that used in j47J to solve the Lindblad equation for the XXZ chain. 
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V. CALCULATING THE CUMULANTS OF THE CURRENT (A QUICK OVERVIEW): 

The principal interest of the Ansatz we constructed is that it allows to calculate the expected values of some 



observables without having to diagonalise explicitly. The question is, then, how we can use formula (27) to obtain 
explicit expressions of the cumulants of the current. As one can see in we did manage to make that calculation. 

However, a certain amount of guesswork was used, and we believe that there is a simpler and more compact way to do 
it than the one we used. For that reason, we will not expose here the full detail of those (rather tedious) calculations, 
but rather a quick overview of the principles behind it. We hope to be able to do the full calculation in an elegant 
and self-sufficient way in the near future, and possibly to do the same for observables other than the cumulants of the 
current (as we said earlier, we should have access to any spatial observable in the s-ensemble). 



A. Expressing of iJ(/i) and /i as parametric infinite series: 



The main point of our reasoning from here on is that we expect the solution to have the same structure as in the 
periodic case |35j . i.e. we expect to be able to write E{ii) and fj, as two infinite logarithmic series in a parameter B 
which goes to with /x: 



oo 



(47) 
(48) 



k=l 



From calculations using our Ansatz in the periodic TASEP, and comparing them to the results from the Bethc Ansatz 
[35j , we were able to determine that this parameter B is proportional to ^^^—^ — - , where we recall that A^ is the largest 
eigenvalue of (J7^T^) (which is the one associated with and goes to 1 when fj, goes to 0) 

Luckily, there is a w 
for. Let us first define 

' ^ 'U,T,-\\P^)){{1\\ 



Luckily, there is a way to write \\Pfj}) as a series in ^^—^ — ^, from which we could get the formulae we are looking 



(l-c-A") 



which is finite for /i — > 0. We can now write (by taking formally the limit A; — > oo in (25)) 



IIP,,)) = 



U.,T, 



(||P*))((l|| + (l-e-A')^) 



A 



A„ 



or equivalently 



K\\P,l 



(1 - ^-n 

A„, 



||P^))-E(^)'I^*) 



\\P, 



(1 



A,, 



which is a well defined series in '^^^—^ — ^ . 



From this, we get 



E{t,)^il\\M^\\P,))=Y, 



{l\\M^{U^T^nP*)) / (1-e-^) 



fc+i 



fc=0 



(l-e-^^) 



A„ 



where iinM^iu^T^np*)) j^^^^^ ^ ^ q because ((1||M^ 



( 

We also get, tautologically (since ((1||P^)) — 1), 

A* = - logfl - (1 - c-^)((l||P^)))l = - log 



which we can then expand in '^^—^ — ^. 



/i. 



i-5^((i||(cv5;)1|P*) 



(l-c- 



fe+i- 



k=0 



(49) 



(50) 



(51) 



(52) 



(53) 
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B. Inferring the final formulae: 



There is a major difference between expressions (53 1 and (48) (or between (52) and (|47|): the coefficients Ck and 

^^lllA-/ (U~^t' )^\\P*')) 

Dk should not depend on /i, but — — " — for instance, does. From this point on, the calculations become 
less precise. The reasoning is as follows: 

• we postulate that the coefficients Ck and should be a somewhat identifiable part of, respectively, ((l||(C/pT^)'^||P*) 
and {{l\\M^{U,,T^nP*))- 

• we then calculate the equivalent terms for the periodic TASEP (which is the simplest case to which our Ansatz 
applies), namely {{lN\\{TP''''f\\lN)) and ((lAr||MP^''(rP^'')'=||ljv)), for A: = 2, using the matrix product formalism 
and the q-deformed oscillator algebra (which, for g 0, becomes a random walk on N); 

• we compare the result with the coefficients from [35], and identify in which part of the calculation they emerge; 



• we isolate the corresponding part of ((l||([/^r^)'=||P*)) and ((1||A/^(C/^T^)'=||P*)), and inject it in ([48]) and ( [47| ; 

• we check numerically (i.e. do exact calculations on small sizes and orders of /i) that our conjecture is correct. 

The results of this calculation for the open TASEP can be found in . The generalisation to the open ASEP [33] 
comes from applying the same reasoning to the periodic ASEP [35]; in that case, we also checked our results against 
DMRG calculations for low order cumulants {E2 to E4) and larger sizes of the system (up to 100 sites). 



VI. CONCLUSION: 



In the present paper, we define and expose the algebraic proof of a matrix Ansatz which gives access to the principal 
eigenvalue and eigenvectors of the current-counting Markov matrix of the open ASEP, for any size and any value of the 
parameters. Using this Ansatz, the author and collaborators were able to obtain exact expressions for the cumulants 
of the current in the open TASEP [32] and the open ASEP [33], which had been an open question for many years in 
the field of non-equilibrium statistical physics. 

Much remains to be done on this subject, and we believe that this Ansatz still has a lot to offer. For instance, as 
was argued in section |III[ the principal eigenvectors of hold the probabilities of observing a given configuration 
conditioned on the current flowing through the system; in other words, it should allow us to analyse the best profiles 
(in the sense of most probable) to produce a given atypical current. The question of finding the optimal path to 
produce a rare event is an important one, notably in the context of complex chemical reactions, and much work has 
been done to find algorithms that produce this optimal path [33] • Obtaining an exact analytical result on that type 
of problem could provide valuable insight or help devise simpler algorithms. 

Another situation where our method could be of use is the Symmetric Exclusion Process, for which many results 
are known for the large size limit, and have been obtained using a coarse-grained description of the system named 
'Macroscopic Fluctuation Theory' [iS] and the related 'additivity principle' [7] [5], but only a few exist for finite sizes 
[46] , The limit g — J> 1 cannot be taken directly in our results, but the present Ansatz can still be applied to the 
symmetric case with a few crucial alterations. However, we have yet to analyse that case in detail, which we intend 
to do in the near future. 

There remains also the question of determining how specific or general the method we have applied here could be. 
We have shown that by defining the transfer matrix U^T^j, which commutes with M^, and then taking jj, to 0, one can 
retrieve the original matrix Ansatz |33| . It would be interesting to know whether the same procedure can be applied 
to other models with matrix product eigenstates, and if there is anything general or generalisable to it. We would 
also like to have a clear idea of the physical significance of that transfer matrix, which we used as a mere calculation 
tool, but might be interesting in itself. 

The author would like to thank K. Mallick and R. Vasseur for useful discussions and for reading the manuscript, 
as well as F. van Wijland and V. Pasquier for helpful comments. 
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Appendix A: Gallavotti-Cohen symmetry for the currents 



In this section, we prove that two current-counting Markov matrices Afj^.j and M^^'y ar similar (and therefore 

have the same eigenvalues) as long as X]^=o l^i ~ ^f=o l^i- 

Let us consider the diagonal matrix Ri{vi) (with 1 < i < L) which multiplies by e"^' all configurations for which site 
i is occupied. The transformation Ri{vi)~^Mi^^.yRi{vi) acts only on Mi_i(/ii_i) and Mi{p,i), and a trivial calculation 
gives 



Ri{vi) ^Mi_i{n,_i)R,{vi) = Mj_i(^i_i - V,) 



(Al) 



and 



R,{viy^Mi{^ii)R,{vi) = M,{fi, + Vi) 



(A2) 



One can therefore transfer any fraction v oi ^ ^ X^iLo l^i from one bond to the previous or the next one. Using 
this, one can go from {/i^} to {/ij} step by step, or simply derive the global similarity matrix such that 



which one can easily find to be 



(A3) 



■1=1 j=0 

There is a particular set of weights {/ii} defined by 



(A4) 



{no = A log 



= A log {q), flL = A log 



(A5) 



for which M^^.j becomes 



A^o(mo) 



-a -y^^^a^ 





-q q^ 

q^-^ -1 





MM 



(A6) 



which is the deformed Markov matrix measuring the entropy production. We sec immediately that 



(A7) 



which proves the Gallavotti-Cohen symmetry for the eigenvalues and between the left and right eigenvectors of AI\ 
with respect to the transformation X ^ (1 — A). 



Considering that fi — A log ( ^ 
namely 



we also obtain the Gallavotti-Cohen symmetry related to the current. 



E{fi) = E[ log 



ap 



(A8) 



which is also valid for the other eigenvalues of Af^, and the corresponding relations between the right and left 
eigenvectors. 

The equilibrium case, where — 1, is somewhat pathological: there is no entropy production, and all the odd 

cumulants of the current vanish. However, by expanding all these relations to second order in the small parameter 

log (S,-'), o„ ca„ ™ .he fl„c.u..io„.di.ip..i„n 

For more results on the Gallavotti-Cohen relation, one may refer to |12) . 
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Appendix B: Derivation of equation (21 1 for the open ASEP 



The derivation of the commutation relation (21) can be carried out in two steps: first, we express the commutator 
of Mi (for 1 <i < L) with either or (the two results are similar) using (28), and show that for both, the sum of 



those commutators cancels out except for two terms, related to each of the boundaries. Secondly, we check that those 
boundary terms, as they appear in the commutator of with the product J7^T^, cancel out as well, using (29) and 



For convenience, we will here write [/„ and r„ as: 



^-{W\A,f{x^^^\V) 



Zl 



i=l 



T,^{W\A,\{X^V). 



(Bl) 
(B2) 



with 



X 



(0 _ 



1 e 
d 1 



(B3) 



Let us therefore consider [Mi, J/^] and [M^, T^]. The elementary matrix Mi acts only on sites i and i + 1, so that 
we only have to consider the commutation with the part of the matrix products that correspond s to those sites. In 
both C/p and T^, that part is Let us write its components in the same basis as Mi in (14), and recall the 

expression of Mi. 



■ 1 e e ee 

d I ed e 

d de 1 e 

dd d d 1 



M,, 



0" 

-(? 1 
0^-10 





(B4) 



We now calculate: 



• 

(1 — q)d de — q 1 — q ed (1 — q)e 

{q — l)d q — de q ed ~ 1 (q ~ l)e 

. 

• 

{1 — q)d de — q ed {1 — q)ed {1 — q)e 
(q — l)d {q — l)de q ed — de {q — l)e 




(1-9) 



0- 

q{ed - 1) 1-ed 

q{l -de:} de-1 















" 








■ 


d 


1 


ed 


e 


-d 


—de 


-1 


— e 


. 












(B5) 



by using the first relation in (|28|) to go from the second line to the third. 

(1-9) 



Let us now define the 'hat' matrices X by X^^t' = (— 1)'^ Xt.t' 



(1-9) 


■ 1 e ■ 


_ (1-9) 


'1 " 


2 


-d -1 


2 


-1 



x 



(B6) 



(where we note . the product in the 2-dimensional space corresponding to the occupation number on one site). 
If we replace the first or the second X in X'^^X'*"*"^^ by X, we get: 



(1-9) 
2 



1 e e ee' 

d 1 ed e 

—d —de —1 —e 

~dd —d —d —1 



(1-9) 
2 



1 e e ee' 

—d —1 —ed — e 

d de 1 e 

—dd —d —d — 1 



(B7) 
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and we finally find that 



(B8) 



The relation equivalent to this one in the matrix product formalism (i.e. after making multiple tensor products) 
can be related to the one found in j37j and used to derive the matrix Ansatz for the multispecies periodic ASEP. 

Putting this relation back in [/^ or and summing over i cancels out all the terms except for those containing X^^^ 
and X^^\ because all the other X*^*' appear exactly twice (once in [Mi, Up] and once in [Mi+i, C/^]) with opposite 
signs. Ultimately, we get: 



L-l 
i=l 



1 L L-l 

1=2 7=1 



L-l 

i=l 



L-l 



1=2 



i=l 



which we may write as: 



L-l 



i=l 
L-l 



i=l 



SO that: 



L-l 



— fl^^'>T — fj^'^'^T + JJ T^^'> — JJ T^-^) 



(B9) 
(BIO) 

(Bll) 
(B12) 

(B13) 



We now need to check that [Moi^i), U^^T^] = -Ul^^Tp ~ Uf,T^'' and [Ml, U^,T^] = Ul^^'T^ + U^.T^''' . As before, Mq 
acts only on site 1, so that only X^^^^ is affected, and the same goes for Ml and site L. 
Let us recall: 



=1(1) 



^L) 



Mo(m) 



—a 7e ^ 
ae'' —7 



Mr 



-s p ■ 

5 -/3 



(B14) 



We calculate: 



[Mo(m),X(i)] 



7e ^d — ae^e (7 — a)e 
(a — 7)d ae''e — 'ye^^^d 



/3d -(5e (/3-(5)e' 
((5-/3)d Se-pd 



(B15) 



By projecting these equations on, respectively, and [V), we get: 

Mo{^l)Aw[ApX^'^] - {w[ 



(70? -ae)^^ A^(7-a)e 
A^(Q!-7)d {ae-^d)Af^ 

I3d-Se il3-d)e 
{5-P)d 5e~pd 



\V) 



(B16) 



(where we used the second and third relations in (28) to get rid of the /is). We naturally find the same expressions 
for (1^1 and [V). 



We can then use relations (|29| and (|30|) to simplify those four equations. We obtain: 

{\-q){W[Ap 





= (a 


-7)(W^I^M 


Mo(Ai), (VK|A,,X(1)' 


= 


-l){W\Ap 






= (/3 


-S) 


■-1 e" 
~d 1 






= (/3 


-S) 


■ 1 e 

-d -] 



1 -e 
d -1 



-1 
1 



-1 -e 
d 1 



\V) 



1 

-1 



1^) 



(B17) 
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so that, on site 1: 



Mi 



+ ((W^|A^X(i)).((a-7)(iy|A 



-1 -e 
d 1 



(B18) 



Seeing that 



1 -e 
d -1 



1 

-1 



and that 



this equation cancel out. 

What's more, we can write (1 — q) 



-1 -e 
d 1 



-1 
1 



X, the first and third part of the right side of 



"-1 0' 


_ (1-9) 


■-1 0" 


1 


2 


1 



X+^-^X. 



"1 
1 



so that 



(1-9) 



-1 
1 



(B19) 



which means precisely that [Mo{ii), Uf^T^] = —Uli — C/pT)i . The exact same calculations on the other boundary 
lead to [Ml, C/^T^] = ul^^T^, + U^,fl^\ and this concludes the proof. 



Appendix C: Derivation of equation (21) for the periodic case 



The periodic case is much simpler than the open one. Equation (21 ) takes the form: 



(CI) 



with 



i=l 



(C2) 



By the same calculations as in the previous section, we find: 



L-l 

1=1 



(C3) 



and we need to check that [Ml{h),X^^^'^ A^X^^^\ = X^^^A^X^^^ ~ X^^^ A^X^^^ 
We have: 



A^d 



Af,e 



eAfj,e 



Af^ eAfj^d &A^ 



dAi_^ dA^G -^f^ 
dA^d dA^ A^,d 



Af,e 
A,. 



0- 

-q e^ 

qc-f" -1 





(C4) 
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so that: 



e^'dA, 



-^9 Af,d - dA^, 















dAfj,e e 





e^Ap - g eAf,d 








q eAf, 



A^e 



q{e-^'eA|, - A^e) 

qie-f'eA^d- A^) e"^^ 

g(e-^A^ - dA^e) e^dA^e 

q{e-^'Af,d-dA^) e^dA^ - 



-g)A,,d A^(de-q ed) {1 - q)eA^d (1 
-l)dA^ (q-l)dA^e Ap(ged-de) (q 



-e^^ 0- 

eA^d 

~A,, 

A^d 0. 











(1-9) 









A,. 



-dA^ -dAf^e 





A, 
-A,. 



eA^d 





■ 

eA^ 
-A^fi 





q)eA^ 
^)A,t 




(C5) 



(we used the second and third relation from ( 28 ) to get from hne one to hne two, and the first to get from line two 
to line three). 

Let us note that this same calculation can in fact be used anywhere in the bulk of the open system in order to 
validate the Ansatz using equs. (34). 



Appendix D: spin-| XXZ chain with nondiagonal boundaries 

In this section, we explain how our construction for the open ASEP can be translated for the spin-^ XXZ chain 
with non-diagonal boundary conditions |49j and hint at a possible relation with the recent solution of the XXZ chain 
with a Lindblad boundary drive [4 7) . 

Let us first define the bulk Hamiltonian of the XXZ spin chain of length L: 



L-l 



(Dl) 



k=l 



with hi acting as: 



hi = 



AO 0" 

-A 1 

1 -A 

A 



(D2) 



on sites i and i + 1 (in the same basis as for cqu. (jll))), and as the identity on the rest of the chain. 

Let us then consider the von Neumann equation for the density operator p with the XXZ Hamiltonian with boundary 
terms h^ and h^ acting only on sites and L: 



ih 



dp 

dt 



[H, p] with H = ho + Ht + hL 



(D3) 



dp 

dt 



We will now show how our matrix product construction can be used to find a solution of the stationary equation 
0. 



Let us write the deformed Markov matrix M^j^^.y for the special choice of weights defined by: 



{fj.0 = - log 



+ t9o, pi^^\og{q), HL^^log 



(D4) 
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which is on the Hne /i = 5 log ( ^ ) +iM. and for which M^.y is Hermitian. The deformed local matrices become: 



—a y/a^ e 



/a-f e 



0- 

-q ^ 





MM 



(D5) 



It is straightforward to check that in this case, we have Mj^;.} = + e, where e is a constant, with the boundary 
matrices being equal to: 



(1- 17-0; + 7) 2^/076"**" 
2^076*^0 (-1 + 9 + a -7) 



2^9 



(-l + g + /3-(5) 2y/j5e^^ 
2v^e-*''^ {I -q- 13 + 5) 



(D6) 



The transfer matrix U^^^T^^^ defined in (34) is therefore a solution to equ. (D3). It might not be a suitable density 
matrix, as its eigenvalues might not be positive, but we may in any case define one by: 



UT{UTy 



Tr 



UT{UTy 



(D7) 



(where we omitted to write the dependence in 
We can also rewrite expressions ( 34 1 in a way better suited to this situation: 



U{C,C') 
T{C,C') 



1 



i=l 

L 



with 



and 



Y = 



1=1 



TV 



N 



N = A„ 

2 2 

S_ ^ Ati±dAfH 
2 2 

(01 - (W^l^(^„-i^) 



(D8) 
(D9) 

(DIO) 

(Dll) 
(D12) 
(D13) 
(D14) 
(D15) 
(D16) 
(D17) 



Matrices iV, 5*+ and S- satisfy a special parametrisation of the Uq[SU{2)] algebra [5T, which is the bulk symmetry 
of the XXZ chain: 



S+N = —NS+ 

NS- = —S-N 
%/9 



(D18) 
(D19) 
(D20) 



It was surprising to find that this solution has a structure almost identical to that of the Lindblad master equation 
found in [47] . where Y is noted and Y is noted S. In that case, it seems that p can be written in the simpler form 
p = y^^^t] ; where the algebraic relations satisfied by the boundary vectors (01 and \ip) might be different from ours. 
It would be interesting to understand the precise relation between those two a priori very different situations. 
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